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Abstract 



A search for a dark matter (DM) annihilation signal into 7-rays toward the 
direction of the Canis Major (CMa) overdensity is presented. The nature of CMa 
is still controversial and one scenario represents it as a dwarf galaxy, making it 
an interesting candidate for DM annihilation searches. A total of 9.6 hours of 
high quality data were collected with the H.E.S.S. array of Imaging Atmospheric 
Cherenkov Telescopes (lACTs) and no evidence for a very high energy 7-ray sig- 
nal is found. Upper hmits on the CMa dwarf galaxy mass of the order of 10^ 
M0 are derived at the 95% C.L. assuming neutrahno masses in the range 500 
GeV - 10 TeV and relatively large annihilation cross-sections. Constraints on 
the velocity- weighted annihilation cross section {(jv). are calculated for specific 
WIMP scenarios, using a NFW model for the DM halo profile and taking ad- 
vantage of numerical simulations of hierarchical structure formation. 95% C.L. 
exclusion limits of the order of 5 x 10~^^ cm^ s~^ are reached in the 500 GeV - 
10 TeV DM particle mass interval, assuming a total halo mass of 3 x 10® Mq. 

Subject headings: gamma rays - dark matter - galaxies: dwarf 



1. Introduction 

It is widely believed that around one third of the energy content of the Universe is 
made of cold matter, of which only a small fraction is luminous and is of baryonic origin. 
In the framework of the Cold Dark Matter (CDM) scenarios, most of the matter is com- 
posed of non-baryonic Weakly Interacting Massive Particles ( WIMPs) . The standard model 
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of particle physics does not provide a natural and suitable candidate for the Dark Matter 
(DM) particle and many theories beyond the standard model hav e been proposed to ex- 
plain its origin and properties (see teertone. Hooper and Silk 2005 ) for a recent review). In 
various models, the self-annihilation of WIMPs gives a 7-ray continuum emission resulting 
from the hadronization of primary annihilation products. The spectral features of such a 
DM annihilation radiati o n might help to distingu ish it from ordinary astrophysical sources 
( iBergstrom et al. 2005al : iBergstrom et al. 2005bl ). The indirect detection of DM through 
very high energy (VHE) 7-rays is one of the best ways to probe the astrophysical nature 
of the DM and the H.E.S.S. array of Cherenkov telescopes, dedicated for detection of VHE 
7-rays in the 100 GeV-10 TeV energy regime, is an interesting instrument for this purpose. 
Many astrophysical objects, ranging from DM clumps to galaxy clusters are expected to 
lead to DM particle annihilation signals detectable with sufficiently sensitive instruments. 
Regions of high concentration of DM are good candidates to search for such annihilations 
and the Galactic Centr e (GC) was first considered. H.E.S.S. observations of the GC region 
( lAharonian et al. 20041 ) revealed a source of VHE 7-ray emission (HE SS J1745-290) but ruled 
out the bulk of the signal as of DM origin (lAharonian et al. 2006d ). Prospects for indirect 
detection in the elliptical galaxy M87 at the center of the Virgo clust er were also investigated . 
The time variability of the VHE 7-ray signal observed by H.E.S.S. (lAharonian et al. 2006al ) 
gave clear evidence that the signal was not of a sole DM origin. There are also other 
candidates with high DM density in relative proximity that might lead to detectable DM 
annihilation signals. Satellite dwarf galaxies of the Milky Way (MW) such as Sagittarius, 
Draco o r Canis Major are popular targe ts, owing to their relatively low astrophysical back- 
ground (lEvans. Ferrer and Sarkar 20041 ). Indeed, dwarf spheroidal galaxies usually consist 
of stellar populations with no hot or warm gas and no cosmic rays, and are among the most 
extreme DM-dominated environments. A null result concerning the search for DM toward 
the Sagittariu s dwarf spheroidal gala xy (Sgr dSph) direction was published by the H.E.S.S. 
collaboration (lAharonian et al. 20081 ). Constraints on the parameter spaces of two popular 
WIMPs models, namely the R-parity conserving Minimal Supersymmetric extension of the 
Standard Model (MSSM) and Kaluza-Klein (KK) scenarios with K-parit y conservation, were 
derived. A null result was also est ablished by the MA GIC collaboration (lAlbert et al. 20071 ) 
and the WHIPPLE collaboration (jWood et al. 20071 ). when searching for a DM annihilation 
signal toward the Draco dwarf galaxy. Upper limits on the velocity-weighted annihilation 
cross-section of DM particle were derived in the framework of minimal SUper GRAvity 
(mSUGRA) models and MSSM models, respectively. 

The present paper reports the search for a DM annihilation signal towards the direction of the 
CMa overdensity with the H.E.S.S. array of Cherenkov telescopes. The paper is organized as 
follows: in Sec 2 the controversial nature of the CMa overdensity is briefiy discussed; in Sec 3 
the analysis of the data is presented, while in Sec 4 the predictions for DM annihilation into 
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7-rays in the CMa overdensity are discussed. Constraints on the WIMP velocity-weighted 
annihilation rate, as well as on the CMa total mass, are given. 



2. A Galactic warp or the relic of a dwarf galaxy? 



Since its discovery (IMartin et al. 2004 ). the nature of the Canis Major (CMa) overden- 
sity is the subject of many discussions o yer whether it is a dw arf galaxy or simply a part of 
the warped Galactic disk. According to (IMomany et al. 20061 ). the CMa overdensity simply 
reflects the warp and flare of the outer disk, a structure frequently observed in spiral galax- 
ies such as the MW. The comparison of the kinematics of the CMa stars with those of the 
Galactic thick disk shows that the CMa stars do not have peculiar proper motions, which 
cast some doubts on the dwarf Spheroidal (dSph) nature of the CMa overdensity. The second 
scenario, which is of interest for the aim of this paper, considers this elliptical overdensity 
as the re mnant of a disrupted dwarf galaxy that could have created the Monoceros "ring" 
structure (IMartin et al. 2004 ). Indeed, numerical simulations show that such a structure 
can be explained by an in-plane accretion event, in which the remnant of the dwarf galaxy 
would have an orbital plane close to the Galactic plane. Another argument in favor of this 
scenario is that the CMa star populations do no t exhibit the same pr operties as those in the 
disk since they have a relatively low metallicity (jde Jong et al. 2007 ). In this case, the CMa 
star population woul d not belong to the Galactic disk, in contrast to what was found in 
(IMomany et al. 20061 ). The mass, luminosity and characteristic dimensions of CMa appear 
quite similar to those of the Sgr dwarf galaxy. As for many dSph, the CMa overdensity 
would thus be an interesting candidate for DM detection. In the remainder of the paper, 
the Canis Major object is assumed to be a dwarf galaxy. 

The CMa over density is located tow ards the Galactic anti-centre direction at roughly 8 kpc 
from the sun (iBellazinni et al. 2004 ) and is the closest observed dwarf galaxy. It is a very 
extended object (Al = 12°, Ab = 10°) with a roundish core approxin iately centered at 1 
= 240° and b = -8° according to various star surveys in this region (IMartin et al. 2004 : 



Martinez et al. 2004 ). In contrast to the Sgr dSph, neither dispersion velocity measure- 
ments, nor luminosity profiles are available so that an accurate modelling of the CMa DM 
halo profile is not possible. However, there are enough constraints to estimate the expected 
7-ray flux from DM particle annihilations in this object. The annihilation cross-section is 
given by the particle physics model (see section IV). As concerns the mass content of CMa, 
the narrow disp ersion between the average mass values found for different dSph galaxies in 



the local group (IMateo 19981 : IWa. 



ker et al. 20071) is an indication that dSph's may possibly 



have a universal host halo mass (IDekel and Silk 19861 ). The mass of the CMa dwarf galaxy 



can then be inferred to be in the same range as the Sgr dwarf galaxy and many other dSph's 
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so that the CMa tot al mass would range between 1 0^ and 10^ flMartin et al. 20041 ) . For 
instance, reference (lEvans. Ferrer and Sarkar 20041 ) gives a model where the CMa mass is 
taken as 3 x 10* Mq. The H.E.S.S. large FoV covers a large part of the CMa core, optimizing 
the chances to see a potential DM annihilation signal. 



3. H.E.S.S. observations and analysis 
3.1. The H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes 



H.E.S.S. is an array of four Imaging Atmospheric Cherenkov Telecopes (lACT's) ( iHofmann et al. 20031 ) 
located in the Khomas Highland of Namibia at an altitude of 1800 m above sea level. The in- 
strument uses the atmosphere as a calorimeter and images electromagnetic showers induced 
by TeV 7-rays. Each telescope collects the Cherenkov light radiated by particle cascades in 
the air showers using a large mirror area of 107 m^ and a camera of 960 photomultiplier tubes 
(PMT's). The four telescopes are placed in a square formation with a side length of 120 m. 
This configuration allows for a accurate reconstruction of the direction and energy of the 
7-rays using the stereoscopic technique. The cameras cover a total field of view of 5° in diam- 
eter. The energy threshold of the H.E.S.S. instrument is approximately 100 GeV at zenith 
and its sensitivity allows to detect fluxes larger than 2 x 10~^^ c m~^ s~^ above 1 TeV in 25 



hours. More details on the H.E.S.S. experiment can be found in (lAharonian et al. 2006bl ). 



3.2. Data processing 



Observations of the CMa dwarf galaxy with H.E.S.S. were carried out in November 2006 
with pointing angles close to the zenith and extending up to 20°. The nominal pointing direc- 
tion was 1 = 240.15° and b = -8.07° in Galactic coordinates. The data were taken in "wobble 
mode" with the telescope pointing ty pically shifted by ±0.7° from the nominal target po- 
sition (IBerge. Funk and Hinton 20071 ). The wobble observing mode reduces the systematic 
effects in the background estimates. The dataset used for image analysis was selected using 
the standard quality criteria, excluding runs taken under bad or variable weather conditions. 
The CMa dataset amounts to 9.6 hours of live time after quality selection. 
Two different techniques are combined for the data processing. The flrst technique computes 
the "Hillas geometrical moments" of the shower images to reco nstruct shower geometr y and 
energy, and to discriminate between 7-ray and hadronic events (lAharonian et al. 20051 ). The 
second technique uses a semi- analytical model of air showers which predicts the expected 
intensity in each camera pixel (Ide Naurois et al. 20031 ). Here, the shower direction, the im- 
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pact point and the primary particle energy are derived with a hkehhood fit of the shower 
model to match the images. Standard cuts for 7-ray/hadron separations are derived by 
simulations. Both analyses provide an energy resolution of 15% and an angular resolution 
better than 0.1°. The combination of these two techniques, referred hereafter as "Com- 
bined Hillas/Model analysis", uses a combined estimator (the so-called "Combined cut") 
and provides an improved background rejection. The background is estimated following the 
template background method. This method uses events that fail the 7-ray selection cuts. 
The template background modelling is well suited for the detection of sources positioned 
anywhere in the FoV, as it estimates the backgrou nd in each sk y direction. More details on 



this background subtraction method are given in (iRowell 20031 ). A charge cut of 60 photo- 
electrons on the image size is applied, as well as a cut on the primary interaction depth of 
particles to optimize the signal-to-noise ratio. Events that pass the analysis cuts are labelled 
as "7 candidates" and are stored in the so-called 7 candidate map n'^'^'^^'^^^'^ {I , b) . Events 
that do not pass the analysis cuts are defined as "background events" and are stored in the 
so-called background map nbck(^ b). Table [1] shows the different cut values used to select the 
7-ray events. 

The 2.5° X 2.5° excess sky map is obtained by the following equation: 

n™(Z, b) = n;^^''"'^'%l, b) - ail, b) x nbck(/, &), (1) 



where a{l,b) refers to the template normalisation factor as described in (IRowell 20031 ). To 
search for a gamma-ray signal, the raw fine-binned maps are integrated with a 0.1° radius 
around each point to match the H.E.S.S. angular resolution, resulting in new oversampled 
maps of gamma-ray candidates and background events, and a cor responding gamma-ray ex- 



cess map. Using the prescription of Li and Ma (iLi and Ma 19831 ) to derive the significance 
for each point of the oversampled map on the basis of the gamma-ray candidate and back- 
ground counts and the template normalization factor, no significant excess is found at the 
target position or at other points in the field of view (Fig. [T]left panel). The distribution of 
significances for the entire map is shown in the right panel of Fig. [1] and is fully consistent 
with statistical fluctuations of the background signal. As the excess map does not show any 
signal, an upper limit on the number of gamma-ray eve nts for each point in the rn ap can 



be derived using the method of Feldman and Cousins (iFeldman and Cousins 19981 ). The 
uncorrelated 7 candidate and normalized background maps, plotted on a 0.2° x 0.2° grid to 
have bins not smaller than the H.E.S.S. angular resolution, are used for the upper limits 
calculations. 
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4. Predictions for Dark Matter annihilations in the Canis Major overdensity 

The DM particles are expected to annihiliate into a continuum of 7-rays through var- 
ious processes such as the hadronization of quark final states, hadronic decay of r leptons 
and subsequent decay of mesons. Two DM candidates are commonly discussed in litera- 
ture: the so-called neutralino arising in supersymm etric extensions of the standard model 



(SUSY) (jjungman. Kamionkowski and Griest 19961 ). and the first excitation of the hyper 



char ge gauge boson in Un iversal Extra Dimension theories (UED) called the B^^^ parti- 



cle (jServant and Tait 20031 ). Typical masses for these DM candidates range from 50 GeV 
to several TeV. The value of the annihilation cross-section is constrained to give a ther- 
mal relic abundance o f WIMPs that is in agreement with the WMAP+SDSS derived value 



(ITegmark et al. 20061 ). The velocity- weighted annihilation cross-sections can be as low as 



10 cm^ s ^, for scenarios involving co-annihilations processes, and be as high as 10 

s —1 

cm s . 

The expected fiux cj)^ of 7-rays from WIMP annihilations occuring in a spherical dark halo is 
commonly written as a product of a particle physics term (d$^^/di?^) and an astrophysics 
term {f^^): 



X r (2) 



The expected number of 7-ray is then given by: 

Ar^ = ToNX / A,s{E-,)(l)-,{E^)dE^, (3) 
Jo 

where Tqn denotes the ON-source exposure time, which depends on the pointing direction, 
and Aes{E^) the averaged H.E.S.S. acceptance during data collection. The velocity- weighted 
cross-section for WIMP annihilation {av) and the WIMP mass are fixed to compute the 
particle physics term in Eqj2j 

d<l>PP {av) ( dN 



dE^ 47rm^j^ 



( dE^ ] dm' 



where (dA^/di?^)DM is the 7-ray spectrum originating for DM particle annihilation. The 
shape of the continuum 7-ray spectrum predicted in the framework of the phenomeno- 
logical Minimal Supersymmetric extension of the Standard Model (pMSSM) depends on 
the model in a complicated way. A simplified parametrization of this shape, for higgsino- 
l ike neutralinos mainly annihilating via pairs of W and Z gauge bosons, was taken from 



( iBergstrom. UUio and Buckley 19981 ). In the case of KK B*^^^ particle annihilations, the 
branching ratios to final states are independent of the WIM P mass. The differential p hoton 
continuum has been simulated with the PYTHIA package (jsjostrand T. et al. 2003 ) using 
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branching ratios from (jServant and Tait 20031 ). 



The astrophysics term /^^ illustrates the DM concentration dependency of the expected 
7-ray flux toward the pointed source: 



= / / p^(/)d/dfi, (5) 
J An J\os 

where p{l) is the mass density profile of the CMa dwarf galaxy and Af2 the detection solid 
angle (Af2 = 10"^ sr, corresponding to the integration radius of 0.1°). 



4.1. Model of the Canis Major Dark Matter halo within the ACDM cosmology 

The purpose of this section is to explain how the astrophysical term f^^ (Eq. [21 Eq. [5]) 
was calculated as a function of the total CMa dark halo mass. The estimate of the astrophys- 
ical term /^^ relies on the modelling of the CMa DM mass distribution. Observationally, 
the DM mass content of dSph galaxies can be derived using velocity dispersion measure- 
ments of their stellar population as well as their luminosity profile. The comparison between 
models and observations can constrain the parameters of their assumed density profiles. In 



of its density profile in the same way as in the literature ( 


Evans. Ferrer and Sarkar 2004: 


Colafransesco, Profumo and UUio 2007: 


Aharonian et al. 2008). 



In the absence of observational data, a standard cusped NEW halo (INavarro. Erenk and White 19971 ) 
was assumed to model the CMa dwarf mass distribution: 

where po is the overall normalisation and the scale radius. The parameters po and 
determining the shape of the profile as well as the halo virial mass M^i^. are found by solving 
the following system of 3 equations: 

M,ir= / Pcuspcd(r)dV (7) 
Mvi, = yp2oo X Rl,, (8) 

where i?vir is the halo virial radius. The virial radius is computed given the virial mass 
Mvir and is defined as the radius within which the mean density equals p2oo (p2oo = 200 x pc, 
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where pc is the critical density of the universe 0). 

Eq. [9] relates the concentration parameter Cvir to the virial mass. The concentration param- 
eter is defined as the ratio between the virial radius and the scale radius R^i-c/r^ in the case 
of a NFW profile. The relation between Cvir and M^--^ is not well-known and was studied 
in various simulations of stru cture formation. Eq. [9] is used following the halo concentra- 



tion fit of (IDolag et al. 20041 ) which is in good agr eemen t with most of N-body simulations 



proposed in the literature (see (IBuUock et al. 20011 ) and (lEke. Navarro and Steinmetz 20011 ) 



as examples). In Eq. [HI z denotes the redshift and h the present day normalized Hubble 
expansion rate, cq and a are the parameters of the halo concentration fit and depend on the 
cosmological scenario (cq = 9.6 and a = —0.1 in a ACDM cosmology). 
The CMa dSph galaxy is located close to the Galactic disk and suffers from strong tidal dis- 
ruptions. A reasonably good estimator of its total dark halo mass is then the mass enclosed 
inside its tidal radius rather than its virial mass. The tidal radius of the CMa dwarf galaxy 
is calculated via the Roche criterion: 

Mdsph(n) _ Muwjd - rt) 

where d is the distance of CMa to the center of the MW. MMw(r) denotes the mass of the MW 
galaxy enclosed in a sphere of radius r. A NFW profile for the MW halo is considered with a 
concentration parameter equal to 10 and a virial mass of 10^^ Mq. The total mass of the dSph 
galaxy is computed by iterative tidal stripping, inserting first Mvir in Eq. [TU] and computing 
successively the total halo mass (using Eq. [7]) and the tidal radius (using Eq. [10]), until the 
convergence of the procedure is reached. The question is now whether or not tidal forces 
significantly remodel the internal structure of tidally affected dSph. Discrepant results hav e 



been reported in the literature regarding this question (IReed et al. 20051 : IStoehr et al. 20021 ). 
Here, it is assumed that tidal forces do not affect the inner part of the density profile so that 
the initial halo structural parameters are kept constant during the stripping procedure. The 
remaining mass is typically found to be an order of magnitude lower than the virial mass. 
The astrophysical term f^^ can then be computed as a function of the halo mass by per- 
forming the line-of-sight integration of the CMa dSph squared mass density, according to 
Eq. O Table [2] shows the obtained NFW structural parameters for a sample of three dark 
matter halos with different virial masses. The integral of the squared mass density increases 
with the halo mass and the f^^ value for a dwarf galaxy of mass 10^ Mq located at an 
heliocentric distance of 8 kpc is found to be in the right order of magnitude (/^^ ~10^^ 



GeV^ cm ^) compared to (lEvans. Ferrer and Sarkar 20041 ). 



^ http://pdg.lbl -gov 
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4.2. Halo independent constraints on the annihilation signal 

Using the upper limit 7 map on the number of 7-events derived with the 7 candidate 
map and the normalized background map (see Sec 3.2 for details), upper limits on the value 
of {av) X can be derived (Eq. El Eq. [3l Eq. ^ for various neutralino masses. As an 
example, the left panel of Fig. [2] shows a sky map of the 95% CL upper limit values on {av) 
X /^P, computed within a pMSSM scenario for a 1 TeV higgsino-like neutralino annihilating 
in pairs of W and Z gauge bosons. The sensitivity is better toward the center of the FoV, 
because the H.E.S.S. acceptance decreases on the edges of the camera. H.E.S.S. is sensitive 
to {av) X values of the order of 10^ GeV^ cm~^ s~^. Since the core of the CMa dSph 
can be anywhere in the FoV, one obtains a distribution of the {av) x f^^ upper limits over 
the FoV. As shown by the right panel of Fig. [2], this distribution is well fitted by a Gaussian 
curve. The mean value can be taken as the {av) x f^^ 95% CL upper limit reference value 
and the corresponding la error bar reflects the uncertainties associated to it. 
Using this procedure, an exclusion curve in the plane {{av) x f^^^ ?tj,dm) can be derived. The 
corresponding curve is shown on the left panel of Fig. [31 The grey shaded areas represent 
the la uncertainties associated with the la error bars on the {av) x f^^ 95% upper limit. 



4.3. Sensitivity to the CMa mass 

As the f^^ factor is related to the CMa halo mass (see Sec 4.1), an upper limit on 
the total CMa mass can be obtained using the previously derived exclusion curve on the 
{av) X quantity and assuming a flxed value for the annihilation cross-section. WIMPs 
velocity-weighted annihilation cross-sections are expected to be of the order of weak-scale 
interaction cross-sections. Exclusion curves in the plane (M^^Jf 

■.fn-DM) are then plotted for 

different annihilation cross-sections within pMSSM scenarios, using the parametrization of 



(IBergstrom. UUio and Buckley 19981 ) for the 7-ray annihilation spectrum. The correspond- 
ing curves are shown on the right panel of Fig. [31 The grey shaded areas represent the 
la uncertainties associated with the la error bars on the {av) x /^^ 95% upper limit, as 
described in the previous section. Annihilation cross-sections larger than 10~^^ cm^ s~^ are 
considered here. Lower cross-sections would have been too small to constrain the CMa dSph 
mass, e.g. the excluded masses would have been of the order of a MW-sized galaxy for 
typical velocity-weighted annihilation cross-section of 3x 10^^^ cm^ s~^. 
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4.4. Sensitivity to the annihilation cross-section of WIMP candidates 

In this part, the CMa total mass is fixed to be 3 x 10® Mq and the corresponding value 
of the astrophysical contribution /^^ in the expected 7-ray fiux is computed following the 
procedure described in Sec 4.1 Table [3] compares the value of the CMa dSph and Sgr dSph 
structural parameters § assuming a NFW profile. The contribution of the astrophysical term 
f"^^ is larger for the CMa dSph than for the Sgr dSph because CMa is closer to the sun. 
Limits on the velocity- weighted annihilation cross-section ^^"^^ss^c.l ^^^^^ then be derived as 
a function of the DM particle mass. Limits are computed in the framework of SUSY and 
KK models. 



The SUSY parameters are computed with the micrOMEGAs vl.37 software package (IBelanger et al. 2004 ) 
Phenomenological MSSM scenarios have been considered. They are characterized by 7 in- 
dependent parameters: the higgsino mass parameter fi, the common sfermions scalar mass 
mo, the Higgs fields vaccum expectation value ratio tanjS, the gaugino mass M2, the trilin- 
ear couplings At and and the CP-odd Higgs mass Ma. Table H] summarizes the region 
of the pMMSM parameter space scanned to generate the models. The left panel of Fig. 
m shows the H.E.S.S. exclusion limits on the velocity weighted cross-section. The black 
points illustrate the computed pMSSM scenarios and the red points r epresent those satisfy - 



ing the WMAP+SDSS constraints on the CDM rehc density QcBMh"^ (ITeemark et al. 20061 ). 
^CDM^^ is allowed to range between 0.09 and 0.11. The H.E.S.S. observations of the CMa 
dSph allows to exclude velocity weighted cross-sections of the order of 5 x 10~^^ cm^ s~^, 
comparable with those derived for the Sgr dSph modelled with a cusped NFW profile. The 
limits obtained are an order of magnitude larger than the velocity-weighted annihilation 
cross sections of higgsino-like neutralinos. 

In the case of KK scenarios, predictions for the veloc ity-weighted cross-section are computed 



with the formula given in (IBaltz and Hooper 20051 ). The expression of (erf) is inversely 



proportional to the squared mass of the lightest Kaluza-Klein (LKP) particle, namely the 
B*^^^ particle. Considered KK models that reproduce the CDM relic measured by WMAP 
and SDSS require a LKP mass ranging from 0.7 TeV to 1 TeV. The right panel of Fig. H] 
shows the H.E.S.S. limits obtained within these models. The H.E.S.S. observations do not 
constrain the KK velocity weighted cross-section. 



^The NFW structural parameter for the Sgr dSph galaxy were estimated in 
(jEvans. Ferrer and Sarkar 20041 ) using the velocity dispersion measurements of the Draco dSph galaxy 
whereas those of the CMa dSph are derived using numerical simulation results (see Sec 4.1) 



5. Conclusions 



The CMa ovcrdcnsity is the subject of many debates over whether it is a dwarf galaxy 
or the warp and flare of the Galactic outer disk. Considering the flrst scenario, its relative 
proximity makes it potentially the best region for searches of a DM annihilation signal. 
However, the lack of observational data prevents the precise modelling of its density profile. 
Assuming a NFW profile and a mass content of 3 x 10^ M© within its tidal radius, typical 
of dwarf galaxies, H.E.S.S. is close to exclude a few pMMSM scenarios with higgsino-hke 
neutralinos, but does not reach the necessary sensitivity to test models compatible with the 
WMPA+SDSS constraint on the CDM relic density. In the case of DM made of B^^^ particle 
from KK models with extra dimensions, no constraints are obtained. 
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the construction and operation of H.E.S.S. is gratefully acknowledged, as is the support by 
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gramme of the CNRS, the U.K. Particle Physics and Astronomy Research Council (PPARC), 
the IPNP of the Charles University, the Polish Ministry of Science and Higher Education, the 
South African Department of Science and Technology and National Research Foundation, 
and by the University of Namibia. We appreciate the excellent work of the technical support 
staff in Berlin, Durham, Hamburg, Heidelberg, Palaiseau, Paris, Saclay, and in Namibia in 
the construction and operation of the equipment. 
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Cut name 7-event cut value 

Combined cut < 0.7 

Image charge min. > 60 photo-electrons 

Reconstructed shower depth min.(rad. length) -1 

Reconstructed shower depth max. (rad. length) 4 

Reconstructed nominal distance < 2.5° 

Reconstructed event telescope multiphcity > 2 



Table 1: List of cuts used in the analysis. The shower depth is the reconstructed primary 
interaction depth of the particles. The nominal distance is the angular distance of the image 
barycenter to the center of the camera. 





Po 


rs 


rt 


M(r < rt) 




(Mo) 


(lO^M© kpc-3) 


(kpc) 


(kpc) 


(Mo) 


(lO^^GeV' cm-5) 


10^ 


4.7 


0.04 


0.28 


3.9 X 10^ 


0.24 


10« 


1.3 


0.28 


1.17 


3.1 X 10^ 


2.2 


lO^o 


0.39 


2.08 


4.15 


1.9 X 10^ 


12 



Table 2: Structural parameters of NFW dark matter halos and associated / values derived 
with the procedure described in the text, for three different virial masses. 





Heliocentric distance 


Po 




^AP 




(kpc) 


(lO^Mo kpc-3) 


(kpc) 


(lO^^GeV cm-5) 


CMa 


8 


1.1 


0.55 


5.9 


Sgr 


24 


1.4 


0.62 


2.2 



Table 3: Comparison table of the NFW structural parameters of the CMa and Sgr dSph, 
with their associated f"^^ values, assuming for both a total mass of 3 x 10® Mo. 
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pMSSM parameters 
lOOGeV^ ^30TeV 
lOOGeV^ mo ^ITeV 

1.2^ tanf3 ^60 
lOGeV^ M2 ^50TeV 
-3TeV^ A,b ^3TeV 
50GeV^ Ma ^lOTeV 



Table 4: Region of the pMSSM parameter space scanned to generate the models. A set of 
free parameters in the considered range is associated to a pMSSM model. 



239 



240 



241 

I (deg) 



-4-3-2-10 12 3 4 



Fig. 1. — (left) Significance map corresponding to the exc ess map compute d in the analysis 
(see text), calculated according to the Li & Ma method (ILi and Ma 19831 ). (right) Signifi- 
cance distribution derived from the significance map. The solid line shows the Gaussian fit. 
The mean value is 0.01 ± 0.004 and the corresponding variance is 1.000 ± 0.005. 




Fig. 2. — (left) 95% CL. upper limit map on the value of {av) x f derived for a 1 TeV 
neutralino (see text), (right) Distribution of the 95% C.L. {av) x f^^ upper limits logarithm 
derived with the upper limit map from the left panel. The solid line shows the Gaussian fit. 
The mean value is 1.66 ± 0.01 and the corresponding variance is 0.14 ± 0.01. 
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Fig. 3. — (left) Upper limits at 95% CL on the value of {av) x f^^ as a function of the 
DM particle mass in the framework of pMSSM scenarios, (right) Upper limits at 95% CL 
on the CMa total mass versus the DM particle mass for different annihilation cross-sections 
in pMSSM scenarios. The shaded area represents the error bars issued from the la error on 
the {av) X i^^ distribution Gaussian fits (see text for details). 




Fig. 4. — Upper limits at 95% CL on the velocity weighted cross-section as a function the 
DM particle mass in the case of pMSSM (left panel) and KK (right panel) scenarios, for an 
assumed CMa total mass of 3 x 10^ M©. The shaded area represents the la error bars on 
^^y^95%c.L ^ggg details), (left) The pMSSM models are represented by black points, 

and those giving a CDM relic density in agreement with the measured WMAP+SDSS value 
are illustrated by red points, (right) The KK models are represented by the black dashed 
line, and those verifying the WMAP+SDSS constraint on ^cdm^^ are labelled in red. 
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